	REAL(8) FUNCTION eps1(t,h)
	implicit none
	real(8) t,h,c
	if(t.ge.273.15) then
	 c = 0.14
	else
	 c = 0.17
	end if 
	if (h==0) h=0.1
!	eps1 = c*(h**(1./7.))*exp(350./t)
	eps1 = c*((h/100.)**(1./7.))*exp(350./t)
	RETURN
	END FUNCTION eps1